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Abstract 



We propose a modification to the standard statistical treatment of the de- 
tection cross-section uncertainties in the analysis of solar neutrino data. We 
argue that the uncertainties of the energy-averaged cross sections of the dif- 
ferent neutrino fluxes in the same experiment should be treated as correlated. 
We show that the resulting allowed regions for the neutrino oscillation pa- 
rameters are significantly larger than the ones obtained with uncorrelated 
uncertainties. 
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Typeset using REVTeX 



The discrepancy between solar neutrino data |]-|||] and the predictions of the Standard 
Solar Model (SSM) [|[], known as the "solar neutrino problem", represents one of the evi- 
dences in favor of neutrino oscillations (see ||). The statistical analysis of solar neutrino 
data allow to obtain information on the values of the neutrino oscillation parameters. The 
latest solar neutrino data [|1]-|4]] have been analyzed in the framework of two-neutrino mixing 
three-neutrino mixing || and four-neutrino mixing [[10 . 



One of the crucial ingredients in the analysis of solar neutrino data is the calculation of 
the statistical covariance matrix V of the uncertainties. This matrix determines the \ 2 °f 
the fit through the formula 



X 2 = E (< r) - R n P) ) t V ~ X )nn ~ • (1) 



Here i?j exp ' ) is the event rate measured in the j th experiment and Rj is the corresponding 
theoretical event rate, that depends on the neutrino fluxes predicted by the SSM and on the 
properties of neutrinos that modify the flux along the neutrino propagation. In the case of 
two-neutrino oscillations the theoretical event rate depends on the mass-squared difference 
Am 2 = m\ — m\ and on the mixing angle i? (see || ) . 

The covariance matrix V includes the uncertainties of the experimental rates R^ xp ^ , and 
the uncertainties of the theoretical rates Rj hr ^ , quantified by the experimental and theoret- 
ical covariance matrices \/( exp ) and V^ thr \ respectively. Since experimental and theoretical 
errors are independent, the corresponding uncertainties can be added in quadrature: 

V = V (exp) + V {thl) . (2) 

Here we are concerned with the covariance matrix of theoretical uncertainties, V^ thr \ 
that can be written as the sum of two independent contributions, 

y(tbr) = y(cs) + y(ix) ? (3) 

where is the covariance matrix of the errors due to the uncertainties of the detection 
cross sections and is the covariance matrix of the errors due to the uncertainties of 
the neutrino fluxes. Very convenient expressions for these matrices have been presented in 
Refs. [n,5|. These expressions have been used by many authors for the statistical analysis 



of solar neutrino data in terms of neutrino oscillations |9|,p^,3,Fi |l0t^3H T5 



Here we would like to examine closely the expression given in [fTTIfl for the covariance 
matrix of the detection cross section: 

^ = ^£(0'^^ ( 4 ) 

i 

where i?^- is the event rate in the detector j due to the neutrino flux produced in the i th ther- 
monuclear reaction in the sun (i = pp,pep, Hep, Be, B, N, O, F) and Cy is the corresponding 
energy-averaged cross section. The quantity AlnCV; = ACij/Cij is the relative uncertainty 
of the cross section CV, , whose value is given in Ref. P|. 

The expression (|j) is based on two assumptions: 1) the errors of the detection cross 
sections in different experiment are independent; 2) the errors of energy-averaged cross 
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section of the different neutrino fluxes in the same experiment are uncorrelated. Indeed, in 
Eq. (U) the errors of energy-averaged cross section of the different neutrino fluxes in the same 
experiment, _R-* hr ^A lnC^, are added in quadrature to the corresponding diagonal element 

y( cs ) 
jj ' 

We think that the first assumption is appropriate, but the second one is questionable. 
This is due to the fact that the correlations between the uncertainties of each detection cross 
section at different energies is not known. Moreover, the uncertainties of the energy-averaged 
cross section of the neutrino fluxes that have an overlap of energy ranges are certainly not 
uncorrelated. 

Since the correlations of the errors of energy-averaged cross section of the different neu- 
trino fluxes in the same experiment are not known, the correct attitude consists in adopting 
the most conservative approach, assuming a complete correlation. In this case the errors 
must be added linearly: 



This approach has been also recommended in Ref. ||16|| , that is the standard reference for 
the detection cross section in Gallium experiments. Let us notice that it reflects a rather 
realistic possibility: that in which the errors of the cross section in a given experiment have 
the same sign at all energies. 

We have performed a fit of the total rates measured in solar neutrino experiments in 
terms of oscillations between two active neutrinos [y e — > or v e — > u T ) in order to study 
the change of the allowed regions for the neutrino oscillation parameters Am 2 and $ when 
the expression (|^) for is used in place of Eq. (§) . The results of our fits are presented 
in Figs. We have used the latest measurements of the event rates in the Homestake 
U and Super-Kamiokande |§] experiments, and the weighted average of the rates measured 
in the two Gallium experiments GALLEX 0] and SAGE |!|. The values of these rates are 
given in Table I of Ref. [[F|. 

Our calculation of the theoretical event rates iC* follows the standard method de- 



scribed in several papers for matter-enhanced MSW transitions fI7|-|I9|l and vacuum oscil- 
lations |20|Jl9[| . We calculate the MSW survival probability of z/ e 's in the Sun using the 
standard analytic prescription [ pT] , P^| , pT| ,|5| and the level-crossing probability appropriate for 



an exponential density profile |22| , |T7|1 . We calculate the regeneration in the Earth using a 
two-step model of the Earth density profile, that is known to produce results that do not 
differ appreciably from those obtained with the correct density profile. We have used the 
tables of neutrino fluxes, solar density and radiochemical detector cross sections available in 



Bahcall's web page [23]. We have neglected the matter effects that slightly affect the vacuum 
oscillation solutions of the solar neutrino problem, as discussed in [0>0]- m Figs- HH1 we 
have used as abscissa the parameter tan 2 $ that allows a better view of the regions at large 
mixing angles with respect to the usual parameter sin 2 2$. The parameter tan 2 $ has been 
employed in the past in the framework of three-neutrino mixing pi-KBl and its use in the 



framework of two-neutrino mixing has been recently advocated in JT5[ because it allows to 
explore the possible presence of allowed regions with $ > tt/4. 



3 



The allowed regions in Figs. |T]-§| are calculated through the constraint x 2 < Xmin + 
Ax 2 (a), where Xmin is the global minimum of the x 2 , xL = 0-42 for Am 2 = 5.1 x 10~ 6 eV 2 
and tan 2 1? = 1.6 x 1CT 3 in both cases of uncorrelated (Eq. (||)) and correlated (Eq. (|5|)) cross 
section uncertainties. We consider confidence levels a = 0.90, 0.95, 0.99 (the regions inside 
the solid, short-dashed and long-dashed contours, respectively) for two degrees of freedom, 
which give Ax 2 (0.90) = 4.61, Ax 2 (0.95) = 5.99, and Ax 2 (0.99) = 9.21. 

The linear addition in of the uncertainties of the detection cross section lead to larger 
values of the diagonal elements of the covariance matrix V and to smaller values of the x 2 i n 
Eq. (|I|), with respect to the case of uncorrelated uncertainties (Eq. (|)). Since the variation 
of x 2 i n the space of the neutrino oscillation parameters tan 2 Am 2 near the minimum of 
the x 2 is proportional to V~ l (the contribution of the variation of V~ x is negligible near the 
minimum) , 

fc> = 2g (*-" - (V-W (g^^) + g^Wj , (6) 

the slope of the x 2 is smaller in the case of correlated uncertainties (Eq. (|5|)) than in the case 
of uncorrelated uncertainties (Eq. (0j)). Therefore, we expect to obtain larger allowed regions 
in the case of correlated uncertainties, with respect to those obtained with uncorrelated 
uncertainties. 

Figures [I] and show the MSW allowed regions obtained using the expressions (f|) and 
(H) for V^ cs \ respectively. One can see that all the three MSW allowed regions, usually 
called SMA (for Am 2 ~ 5 x 10~ 6 eV 2 and tan 2 $ ~ 10" 3 ), LMA (for Am 2 ~ 2 x 10" 5 eV 2 
and tan 2 $ ~ 0.3) and LOW (for Am 2 ~ 10~ 7 eV 2 and tan 2 $ ~ 0.5), are larger when 
the expressions (Q) with correlated cross section uncertainties is used. The minimum and 
maximum values of tan 2 $ and Am 2 in the 99% CL allowed regions in Fig. [I] (unc) and 
Fig. ^ (cor) are listed in Table | and the difference between these regions is illustrated in 
Fig. ^. One can see that the SMA-cor region is slightly larger than the SMA-unc region. 
The LMA-cor region is appreciably larger than the LMA-unc region, allowing values of 
Am 2 and tan 2 i9 as high as 3.5 x 10 _4 eV 2 and 0.82, in contrast to the upper limits 2.4 x 
10 _4 eV 2 and 0.67 for the LMA-unc region. The larger effects occur in the LOW region. 
In the LOW-cor region tan 2 $ can be as large as 1, corresponding to maximal mixing, and 
2.4 x 10" 8 eV 2 < Am 2 < 3.5 x 10~ 7 eV 2 , whereas in the LOW-unc region tan 2 $ < 0.82 and 
5.6 x 10~ 8 eV 2 < Am 2 < 2.6 x 10" 7 eV 2 . 

Figures f| and |5] show the vacuum oscillation allowed regions obtained with uncorrelated 
(Eq. (fflj) and correlated (Eq. (|5|)) cross section uncertainties, respectively. One can see that 
the allowed regions obtained with correlated cross section uncertainties are significantly 
larger that the ones obtained with uncorrelated cross section uncertainties. For example, 
the 99% CL allowed region at Am 2 ~ 5 x 10~ 10 eV 2 is considerably larger with correlated 
cross section uncertainties (Fig. |[) than with uncorrelated ones (Fig. f|), and the 90% CL 
allowed region in Fig. |5| at Am 2 ~ 2.7 x 10~ 10 eV 2 does not exist in Fig. ^. 

In conclusion, we have proposed a modification of the standard statistical treatment of 
the detection cross-section uncertainties in the analysis of solar neutrino data. We have 
argued that, because of lack of information, a conservative approach should be used in 
the treatment of the correlations between the energy-averaged cross sections of the different 
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neutrino fluxes in the same experiment, leading to the assumption of correlated errors and to 
a linear addition of the uncertainties. We have shown in Figs. [l|-|5| that the resulting allowed 
regions for the neutrino oscillation parameters are significantly larger than those obtained 
with uncorrelated detection cross sections uncertainties. We think that the appropriate 
treatment of the detection cross section uncertainties proposed here will be rather important 
in the future, when more precise solar neutrino data obtained in the Super-Kamiokande, SNO 
P7J, Borexino [28], GNO [EO] and other experiments will be published. 



ACKNOWLEDGMENTS 

We would like to thank A. Bottino for useful suggestions and encouragement. C.G. would 
like to thank P. Krastev, E. Lisi, H. Murayama and S. Petcov for useful discussions on the 
fit of solar neutrino data. 



5 



REFERENCES 



[1] B. T. Cleveland et al, Astrophys. J. 496, 505 (1998). 

[2] GALLEX, W. Hampel et al, Phys. Lett. B447, 127 (1999). 

[3] SAGE, J. N. Abdurashitov et al, Phys. Rev. C60, 055801 (1999), |arXiv: astro 
ph/ 9907 1151 . 

[4] Y. Suzuki (Super-Kamiokande Coll.), Talk presented at the XIX Int. Sympo- 
sium on Lepton and Photon Interactions at High Energies, Stanford, Aug. 1999 
( |http : / / lp99 . slac . st anfor d. edu| ) . 

[5] J. N. Bahcall, S. Basu, and M. H. Pinsonneault, Phys. Lett. B433, 1 (1998), |arXiv:astro- 
ph/ 9805135j 

[6] S. M. Bilenky, C. Giunti, and W. Grimus, Prog. Part. Nucl. Phys. 43, 1 (1999), 



|arXiv:hep-ph/9812360 



[7] M.C. Gonzalez-Garcia, P.C. de Holanda, C. Pena-Garay and J.W.F. Valle, |arXiv:he~p- 
ph/99U6465. 



[8] J.N. Bahcall, P.I. Krastev and A.Yu. Smirnov, |arXiv:hep-ph/9911248 . 
[9] G. L. Fogli, E. Lisi, D. Montanino, and A. Palazzo, Phys. Rev. D62, 013002 (2000), 
|arXiv:hep-ph/9912231 . 



[10] C. Giunti, M.C. Gonzalez-Garcia and C. Pena-Garay, |arXiv:hep-ph /0001101 
[11] G. L. Fogli and E. Lisi, Astropart. Phys. 3, 185 (1995). 



[12] G.L. Fogli, E. Lisi, D. Montanino and A. Palazzo, |arXiv:hep-ph/0005261 



[13] M.C. Gonzalez-Garcia and C. Pena-Garay, |arXiv:hcp-ph/0002186 
[14] A. Friedland, |arXiv:hep-ph/ 0002063; . 



[15] A. de Gouvea, A. Friedland and H. Murayama, |arXiv:hep-ph / 0002064 



[16] J. N. Bahcall, Phys. Rev. C56, 3391 (1997), |arXiv:hep-ph/9710491 



[17] T. K. Kuo and J. Pantaleone, Rev. Mod. Phys. 61, 937 (1989). 

[18] P. I. Krastev and S. T. Petcov, Phys. Lett. B207, 64 (1988). 

[19] P. I. Krastev and S. T. Petcov, Phys. Lett. B299, 99 (1993). 

[20] P. I. Krastev and S. T. Petcov, Phys. Lett. B285, 85 (1992). 

[21] S. J. Parke, Phys. Rev. Lett. 57, 1275 (1986). 

[22] S. T. Petcov, Phys. Lett. B200, 373 (1988). 

[23] J.N. Bahcall, WWW page: [http://www.sns.ias.edu/~inb7] . 

[24] G. L. Fogli, E. Lisi, and G. Scioscia, Phys. Rev. D52, 5334 (1995), |hep-ph/9506350| . 
[25] G. L. Fogli, E. Lisi, and D. Montanino, Phys. Rev. D54, 2048 (1996), |arXiv:hep- 
| ph/9605273| . 

D. Montanino, and G. Scioscia, Phys. Rev. D56, 4365 (1997) 



G. L. Fogli, E. Lisi, 
|hep-ph/970623q . 
[27] J. Boger et al, |arXiv:nucl-ex/9910016 . 

[28] Borexino, F. von Feilitzsch, Prog. Part. Nucl. Phys. 40, 123 (1998). 

[29] GALLEX and GNO, T. A. Kirsten, Nucl. Phys. Proc. Suppl. 77, 26 (1999). 



6 



TABLES 





tan 2 $ m i n 


tan 2 "!? ma x 


A ™min (eV 2 ) 


Am 2 max (eV 2 ) 


SMA-unc 


3.3 x 10" 4 


2.9 x 10~ 3 


3.8 x 10~ 6 


1.0 x 10~ 5 


SMA-cor 


3.0 x 10" 4 


3.2 x 10~ 3 


3.8 x 10~ 6 


1.0 x 10" 5 


LMA-unc 


0.14 


0.67 


6.2 x 10~ 6 


2.4 x 10" 4 


LMA-cor 


0.14 


0.82 


5.6 x 10~ 6 


3.5 x 10" 4 


LOW-unc 


0.45 


0.82 


5.6 x 10~ 8 


2.6 x 10~ 7 


LOW-cor 


0.41 


1.0 


2.4 x 10~ 8 


3.5 x 10~ 7 



TABLE I. Minimum and maximum values of tan 2 d and Am 2 in the 99% CL allowed regions 
in Fig. H (unc rows) and in Fig. || (cor rows). 
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FIG. 1. Allowed regions for the neutrino oscillation parameters Am 2 and tan 2 ?? in the case 
of MSW resonant transitions in the Sun and uncorrelated detection cross section uncertainties 
(Eq. (||)). The regions inside the solid, short-dashed and long-dashed contours are allowed, respec- 
tively, at 90%, 95% and 99% confidence level. 
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FIG. 2. Allowed regions for Am 2 and tan 2 "d in the case of MSW resonant transitions in the 
Sun and correlated detection cross section uncertainties (Eq. (|5|)). The regions inside the solid, 
short-dashed and long-dashed contours are allowed, respectively, at 90%, 95% and 99% confidence 
level. 
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99% CL 




in the Sun. The regions within the solid and dashed contours have been obtained, respectively, 
with correlated (Eq. (]5])) and uncorrelated (Eq. (^)) detection cross section uncertainties. 
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UNCORRELATED CS ERRORS 




- 1 o 
10 10 
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FIG. 4. Allowed regions for the neutrino oscillation parameters Am 2 and tan 2 $ in the case 
of vacuum oscillations from the Sun to the Earth and uncorrelated detection cross section uncer- 
tainties (Eq. (|j)). The regions inside the solid, short-dashed and long-dashed contours are allowed, 
respectively, at 90%, 95% and 99% confidence level. 
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CORRELATED CS ERRORS 
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FIG. 5. Allowed regions for Am 2 and tan 2 i? in the case of vacuum oscillations from the Sun 
to the Earth and correlated detection cross section uncertainties (Eq. (|5|)). The regions inside 
the solid, short-dashed and long-dashed contours are allowed, respectively, at 90%, 95% and 99% 
confidence level. 
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